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Abstract 

The surface conductivity for conduction electrons with a fixed chirality in a topological insulator 
with impurities scattering is considered. The surface excitations are described by the Weyl Hamil- 
tonian. For a finite chemical potential one projects out the hole band and one obtains a single 
electronic band with a fixed chirality. One obtains a model of spinless electrons which experience a 
half vortex when they return to the origin. As a result the conductivity is equivalent to a spinless 
problem with correlated noise which gives rise to anti-localization. We compute conductivity as a 
function of frequency and compare our results with the Raman shift measurement for Bi2Se^. 
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I. Introduction 



Topological insulators (TI) are time reversal invariant systems which obey Kramers the- 
orem jl-4|. Toplogical insulators are characterized by the surface excitations and are de- 
scribed by the Weyl Hamiltonian. Bi 2 Se^ is such a material which is described by the Weyl 
equation with a single Dirac cone which lies below the chemical potential /i and the bulk 
gap. Due to time reversal symmetry the backscattering is suppressed in agreement with 
scanning tunneling microscopy (STM) experiment and theory [6j. Conductivity results are 
less conclusive, due to the presence of the 3D bulk gap \^ or the insulating gap observed in 
thin layers of TI. 

The charge current for the Weyl equation is identified with the spin half operator. As a 



result the charge current and the spin current are related to each other |13J . In the presence 
of impurities, elastic scattering conserves energy but not quasi-momentum giving rise to a 
finite conductivity. 

Transport properties have been calculated building on the relation with the spin-orbit 
model, which belongs to the symplectic ensemble jz] and therefore quantum interference 
gives rise to anti-localization. These results have been confirmed by js] and [9 12]. 

The purpose of this work is to compute the surface conductivity for the Weyl Hamiltonian 
with a finite chemical potential The surface Hamiltonian for a single Dirac cone and a 
finite chemical potential has been obtained in ref. 13| and can be described by the Weyl 
Hamiltonian 14J with a Fermi velocity v ~ 5 x 10 6 — and a chemical potential /i ~ 0.1 eV. 
The Weyl Hamiltonian is characterized by opposite chirality for electrons and holes. For a 
large chemical potential we can replace the Weyl Hamiltonian by a single electronic band 
with a fixed chirality. The model resembles a spinless model in two dimensions, therefore the 
presence of elastic scattering might give rise to localization. The backscattering potential 
for a single scattering is zero but multi-scattering are allowed. Therefore we might expect 
localization . Consequently, a particle can backscatter after two scatterings allowing a 
decrease in the conductivity. The difference between the one band spinless electrons and 
our case arising from the fact that the conduction electrons are obtained after a projection 
of the chiral spinors of the Weyl spin half fermions. The effect of the projection modifies the 
random scattering matrix elements but preserve the property of spin half spinor introducing 
correlations which are determined by the momentum difference between the incoming and 
outgoing electrons. In order to determine the quantum corrections we compute the quantum 
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return probability P qm given by the interference between a closed path and the time reversed 
path. Due to the projection which preserve the spin half vortex at k = the quantum return 
probability P qm vanishes. Using these results we compute the conductivity using the method 



employed for spinless electrons, 
the maximal crossed corrections 



he impurity scattering affects the charge and gives rise to 



If]]. We find that the charge is enhanced giving rise to an 
enhanced transport life time. Due to the angle dependent scattering the maximal crossed 
diagrams change sign and give rise to anti-localization. This effect result is similar with the 
results obtained by [8] in graphene with a zero chemical potential. 

The plan of this paper is as follows. In Sec. 77 we introduce the Weyl model. We construct 
the quantized particle and anti-particle bands imposing the time reversal symmetry. Due 
;o the singularity at k = the unperturbed spinor must be obey the the Pfaflian properties 
2)]. In Sec. Ill we consider the model for a single conduction band .In Sec. IV we compute 
the quantum return probability for the one band model. In Sec. V we construct the Green's 



functions 



15| for the conduction electrons. We compute the ladder and maximal crossed 



diagram . In Sec. VI we compute the static and frequency dependent conductivity. In Sec. 
VII we discuss our results and propose that a good confirmation of the surface conductivity 



can be obtained from Raman scattering 



19]. In Sec. VIII we present our conclusions. 



II. The Weyl model 



The surface state Hamiltonian for a topological insulator of the Bi 2 Se^ family materials 
is given by the Weyl model [13 1 . The presence of a random potential and an electromagnetic 
field modify the Weyl model in the following way: 

S= J dt J d 2 r^(r,t)(id t - eA (r,t))^(r,t) - H 

H = hv [ dW(f) \a\-id 2 - ~A 2 (r,t)) - ^(-ift - |Ai(r,t)) + U sc (r)] *(r), 



where A(t) 



E(n) jnt 
in e 



(1) 



is the external vector potential, Ao(r,t) is the external scalar po- 



tential and U sc {r) is the random potential controlled by the white noise correlation function 
((U sc (k)U sc (k'))) = D sc 5 2 [k + k'\ . The unperturbed Weyl Hamiltonian in the first quantized 
form is given by h[k] = —a 2 k\ + a 1 k 2 - Here h[k] is time reversal invariant and obeys the 
transformation relation n d~ l h\k] , d = h[—k), where = ia 2 K is the time reversal operator and 
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K is the conjugation operator. The eigenstates form a Kramers degenerate pair with the 
property fi 2 = — 1. At the point k = the eigenvectors \u^- + \k)}, \u^~\k)) are degenerate 
and we need to choose a representation where the eigenvectors are orthogonal to each other. 
In order to regularize the problem we will replace the Hamiltonian h[k] = —a 2 k\ + a x k 2 by: 

h[k] = -a 2 ^^- o l k 2 + o 3 nM(k), n ->■ 0. (2) 

We introduce the function 9[M{k)\ which is one for M(k) > and zero otherwise. We 
solve the eigenvalue equation in the presence of the mass term M{k) . For the region 

— * — * 

9[M{k)\ = 1 find that that the eigenvector \u\{k)) corresponds to the positive eigenvalues 
and eigenvector \u 2 (k)) to the negative eigenvalues . For the region 9[—M(k)] = 1 the 
eigenvector \u 2 (k)) corresponds to positive eigenvalues and \ui(k)) correspond to the negative 
eigenvalues. Using a mass term which preserve time reversal symmetry M(k) = —M(—k), 

— * 

we guarantee that the eigenvectors at the point k = are orthogonal to each other and obey 
Krammers theorem. We will consider the limit 77 — )► and neglect the energy corrections 
but we will consider the topological effects caused by the mass term M(k) which preserve 
the time reversal symmetry ft 2 = — 1. We note that for the wrapping case, the mass 
M(k) = (ki + ik 2 ) 3 + (ki — ik 2 ) 3 = 1k\{k\ — 3k 2 ) provide a natural regularization for the 
eigenvectors at k — . The eigenvectors \ui(k)) and \u 2 (k)) are given in terms of the 
multivalued phase phase x(k), tan(x(re) - 



1 r.T*. 
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«i(*)> = o l*> ® [1, -ie ix{h Y , |«2(fc)> = - |fc> (8) [ie-^ h \ if , (3) 
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The phase x0) obey the property x(~k) = 7r + XW- \ u i(k)) an d \ u 2(k)) form a Kramers 
pair. Due to the time reversal symmetry the two eigenvectors are related through the 
Pfaffian matrix, W h2 (k) = (ui(k)\0\u2(k)) with = Wl 2 (k)tiu 2 {k)) and \u 2 (-k)) = 

W 2 1 (k)'d\ui(k)) . The spinors \ui(k)) and \u 2 (k)) will be used to construct the eigenvector 
\u( + \k)) for particles and \u^~\k)) for antiparticles for the entire Brillouin zone. 

\u(-\k)) = 6[M(k)]\u 2 (k)) + 9[M(-k)]\ Ul (k), 
\u^\k)) = 9[M(k)]\ Ul (k)) + 9[M{-k)\\u 2 {k)). 

(4) 

The eigenvector \u^ + \k)) corresponds to the eigenvalue E p = E + = hv\k\ (for particles) 
and \u^~\k)} corresponds to the eigenvalue E a = E_ = —hv\k\ (for anti-particle). The 



eigenspinors chage the sign under a full rotation of 2n. In two dimensions a full rotation 
of 2-7T is equivalent to two consecutive inversions. We introduce the inversion operator P, 
P| u (+)(jfe)) = |u(+)(-jfe)), P|uH(jfe)) = | u (-)(-jfe)) and find that P 2 = -1: 



p2 e ±i x (fc) = p e ^x(-fc) = p e ^(x(fc)+T) = e fW-fe)+T) = e ^(x(fc)+T+T) = e »T e fxW = _ e fx(*) 
P 2 \u (+) (k)} = -\u( + \k)};P 2 \u ( -\k)} = -\u(~\k)) 

(5) 

This shows that the two eigenspinors \u^(k)) behave like a spin half fermion under rotations. 

We will include the Fermi velocity v to describe the physical energy spectrum. In order 
to incorporate the velocity in the Hamiltonian we multiply the spinor operator by y/v. 
The spinor operator \J/(r) is decomposed into the eigenmodes of the unperturbed Weyl 
Hamiltonian. 

tf(r) = v^^e* f [C(fc)« (+) (fc) + B\-k)u ( -\-k)} = y^^e^tf). (6) 

k k 

We choose that particle operators C(k), C\k) and the anti-particle operators B(k), B\k) 
should obey anti-commutation relations: 

[C(k),C\k')} + = [B{k),B\k')\ + = 5 {2 \k-k\ [C(k),C(k')} + = [B(k),B(k')] + = 0. (7) 

The spinor projection ^ a (r) = {^{f)\o) obeys the modified anti-commutation relations: 

\* a (f),*l(?)] + = v8 a ^ 2 \r,P), [W^M^U = 0, (8) 

For a a chemical potential /z = hvk F > {k F is the Fermi momentum) we introduce the 
notation |/x > for the ground state. The particle operator C(k) and the anti-particle operator 
B(k) annihilate the ground state, C(k)\/i) = for E(k) > jj, > 0. We will consider a situation 
where the chemical potential /x measured from the Dirac point has the typical value of 0.01 
eV with a Fermi velocity v of 5 x 10 5 — . Using the eigenspinors given in eq. (4) we find 
that the Hamiltonian in the present of the scattering potential U sc (k — p) is equivalent to 
two coupled bands. 

H = Y^K hv \k\ ~ n)C ] (k)C(k) + (hv\k\ + n)B\-k)B(-k)], 

k 

H D = ^^vU sc (k-^\v ip > p \k^C\k)C(p^ + V^ a \-k 



+V M (k, -p)C\k)B\-p) + V {a > p \-k,p)B(-k)C(p) 



(9) 



where Hq is the unperturbed Hamiltonian and Ho represent the effect of the impurity 
scattering. The spinors structure gives rise to vertex functions for the coupling to the random 
potential in momentum space. The vertex functions are given in terms of the particle and 
anti-particle matrix elements (p stands for particles and a stands for antiparticles): 

v M (k,p) = («<+>(£) |««(p)>, v M (-k,-p) = ( u (-\-k))\ u (-\-p)), 

V M (k, -p) = {u i+ \k)\u^\-p)}, V M (-k,p) = (u^\-k)\u^\p)). 

(10) 

From Eq. (1) we obtain the external coupling of the electrons to the electromagnetic field 
A{f,t): 

H ex \t) = j d 2 r[J 1 (f,t)A 1 (f,t) + j 2 (r,t)A 2 (r,t)]. (11) 

From Eq. (1) we also find that the currents are given by: Ji(r,t) = (— e)v^^ (f)cr 2 ty(f), 
J 2 (f,t) = (-e)^ t (r)(-^ 1 )^(r)- 



Using the linear response theory 20j for the current-current correlation function in the 
presence of an external vector potential Ai(t) = %e~ in * we compute the interband absorp- 
tion. We find that the optical conductivity in the limit of weak disorder and frequencies 
Q > 2/i is given by the universal relation a(Q > 2/x) = 

For a finite chemical potential \i we integrate the anti-particle band 

#(«) = ^2{hv\k\ +2^)B\-k)B{-k) + ^^vU sc (k-^\v M (-k,-p)B(-k)B\-p) 

k k V 

+V M {k, -p)C\k)B\-p) + V {a ' p) {-k,p)B{-k)C{p) . 

(12) 

We perform the integration in the limit where white noise fluctuations are negligible in 



comparison with the chemical potential \JD SC < pi. Under this condition the anti-particles 
can be integrated out. As a result the particle Hamiltonian is modified by the induced term 
5H m± which can be understood as shift of the chemical potential. 

5H** = -f**f*« f-CHk,u) 



2tt / 2tt / 2tt 



*2. V M { %,fi 3 vU sc (k-p)U sc (p-k + q)v vi^^k-f) 



C{k — q,cu). 
(13) 



Here E(&, o>; a) is the anti-particle self-energy induced by the disorder. When the chemical 
potential \i obeys the condition \(E(k) + E(fc, u; a)) — (E(p) + E(p, w; a)) — 2/i| < 2/i we can 
perform the momentum p integration which allows the replacement of U sc (k — p)U sc (p— k+ q) 
with ((U sc (k—p)U sc (p—k+q))) = D sc 5 2 [q\. Consequently, the effective one-band Hamiltonian 
will have a new chemical potential //—>■// — . Therefore for chemical potential /z which 
satisfies [i > ^f- the one-band approximation is justified. 

III. The single band model 

As a result of the discussion in the previous chapter we conclude that for a finite chemical 
potential the surface of the TI can be described by a single band with a correlated potential. 

tf( c ) = ^2(hv\k\ -n)C\k)C{k) + ^^vU sc {k - p^V M {k,p)C ] (k)C{p) (14) 

k k P 

Where V^ p ' p \k, p) = (u^fylu^^p)) is evaluated with the help of explicit form of eigen- 
spinors given in equation (4). We introduce new operators C\k) and C(k): 

C-t(jfc) = C\k) \e^6[M(k)\ + e^ x( %[-M(£)]" 
C(k) = ^ x ^6[M(k)} + e* x ®e[-M(k)]]c(k) 

(15) 

As a result we can express the Hamiltonian in eq. (14) in the following way: 

/f(«) = ^(hvfii -n)c\k)c{k) + Y / Y / vC\k)r ll (k,p)d(f) (16) 

% k p 

In eq. (16) T\\(k,p) represents the scattering matrix. 

r,|(fc,p) = U ac (k-p)F\\(k,p); f||(fc,p) = cos[i( X (fc) -x(p))] (17) 

The angle xip) ac ts as a half vortex, when the momentum is reversed the angle xip) increases 
by 7T, x(-P) = xiP + n - As a result we obtain: t\\(k,-p) = cos[|(x(£) - x(~p))] = 
cos[i(x(fc) - x(P)) - f ] = sin[|(x(fc) - xip))] 



For the backscattering direction p — —k we find that the single particle scattering 
vanishes f \\(k, —k) = 0. 

Next we consider the coupling of the external vector potential to the conduction band 
only: 

H ext (t) = j ^r^OAi^i) + 72(^0^, *)] = (-«;) / ^ / ^ 
C\k,t)C(p,t)[(u {+ \k)\ - a 2 \u^ + \p))A 1 (-(k - $,t) + ( U ( + \k)\a 1 \u (+ \p))A 2 (-(k - p),t)} 



18) 



where the vertex functions are given by: 

Wx{k,p) = cos[^( X (k)+x(p))h W 2 (k,p) = -sin[i(x(£) + X (p))}- (19) 

Using the definition of the current operator Ji(q,t) = 
(—ev) J j ■^^Wi(k,p)C^(k,t)C(p,t)5 2 [k — p — q\ we obtain from the linear re- 
sponse theory the induced current 5(Ji(q, t\) = J^dtiRi^^t — t 1 )v4 1 (— q, t x ) with the 



correlation function R lt i(q,t — ti) given by |20j: 



RiMt-ti) = *i]([^(g,^i(-£*i)]}- (20) 

We perform the Fourier transform of Ri^(q, t — ti) with respect the frequency Q and find 

in 



that the conductivity is given by a(q = 0, fi) = ' Rl ' 1< '^ = cr(Q). 



IV. Interference effects and the consequence for the conductivity 

The quantum conductivity is controlled by the multiple scattering process of closed paths. 
The interference between a closed path and the time reversed paths allows to define the 
quantum return probability P qm - If we denote by the amplitude for a particle to return 
to the starting point and T,9~ the time reversed amplitude, the quantum return probability 
P qm is given by P qrn = (3? + T3?)(3? + T£F)* (Tis the time reversal operator). For time 
reversal invariant systems the quantum return probability is P qm = 4(3?)(^)* and for 
spin half systems with spin-orbit interactions P qm = [171 ] . Since the classical return 

probability P c \ = 2(3?)(£?)* is smaller for time reversal invariant and larger for spin half 
systems with spin-orbit interactions one conclude that quantum interference gives rise to 
weak localization for the first case and anti-localization for the second case. 
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The scattering amplitude to order n + 1 for a particle with momentum k to be backscat- 
tering to —k is is given by ^(k, —k): 

-k) = r„ 0, -k) + v 1 ^n,:-^ + v v , mmvMf-\ + 

Y E(k)-E(p)+ie ^y(E(k)-E(p) + ie)(E(k)-E(l) + ie) 

T \\(k,Pi) rv-^ rn(pi,^j)...r||(p n _i,p n ) i rn(p n ,-fc) 



r E(k) - E(fc) + X - E(p 2 ) + ie)..(E(k) - E(p n ^) + te)^ E(k) - E(p n ) + *e 



From equation (21) we observe that the first order backscattering term is zero Tj|(A;, —k) = 0. 

— * — * 

The time reversal amplitude T3?(k, —k) is obtained by reversing all the momenta 

—k) = n-t k) = r„ (-k, k) + v rjt^mtlA + 

" y E (-k-E(-p)+ie 

r||(-fc,-p)rn(-p,-zlr||(-z;fc) | ^ ^(-fc,-^) 

Y r " + ^X^"*) " ^H) + <c) " Y E(k) - -E(pi) + ze ' 



[EE 



ru(-pi,-^)...rn(-p n _i,-p n ) i T\\(-p n ,k) 



P2 



( E (-k) - E(-p 2 ) + ie)..(E(-k) - £(-p n _i) + iey E{-k) - E(-p n ) + it 

(22) 



(21) 



Using the fact that the eigenvalues obey E(—p) = E(p) and the scattering vertex satisfies 
the relation rii(— pi,— p 2 ) = r\\(pi,p 2 ) we find that difference between the two amplitudes 
comes from the terms outside the bracket. Using the relation x{~ k) = x(k) + tc we obtain : 

r||(fc,^)r||(p n ,-fc) = -m-k,-p)m-p n ,k) (23) 

As a result \(^(k, — k)+Tj?(k, — k))\ =0 and therefore P qrn = 0. Next we compute compute 
the amplitude for the scattering of a particle with incoming momentum k to outgoing mo- 
mentum -k = Q: P qm (Q) = (^(k,-k + Q)+T^(k,-k + Q))(3 r (k,-k + Q) + T3 r (k,-k + 
Q))* Using the properties of the vertex functions: 

f||(fc,p)f||(p,-fc + (^)f||(-fc, -p)f||(-£fc-3)«- 

\ sin 2 [ X (£) - X (P)\ cos[±( X (k) - X (k - Q))] cob[±( X (P) ~ x(p~ Q))] cos[i( X (£ - Q) - x(p~ Q))] 
^-]sm 2 [( X (k)-x(p))][l- ( " 



(24) 
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we compute the amplitude for small angle backscattering: 

Pgm(Q) = (^(k,-k + Q)+T^(k,-k + Q))(,?(k,-k + Q)+T£r(k,-k + Q)y Follow- 
ing the analysis performed for the backscattering result obtained for scalar waves in time 
reversed systems [isl we compute the backscattering amplitude at small angle (small Q). 
For incoming electrons with momentum k \k\ ~ kp which are backscattered at an angle 9 
we introduce fcpsin[0] = \Q\ and find in terms of the the elastic mean free path / [3] : 

P qm {Q) = {^{k, -k+Q)+T3T(k, -k+Q))(&(&, -k+Q)+T^(k, -k+Q))* « ' ' 



[l + k F lsm[6]) 2 
(25) 



It is important to stress that this result has been obtained for a time reversal invariant 
system in the presence of a half vortex at k = 0. 

V. The Green's function for spinless electrons in the presence of disorder 

In order to compute the conductivity we need to include the interference results estab- 
lished in section IV. This will be done using the method of Green's functions for a random 
potential described in the literature jsQ | . We will use the Green's function G(k, s;p, 0) for 
the conduction band in the presence of the random potential and G^(k,s), the Green's 
function in the absence of disorder: 

G(k,s;p,0) = -i(T(C(k,s)C\p,0)), 
G {0 \k, s) = -t(0\T{C{k, 5)^(^,0)10), 

/oo 
dse icJS G^{k, s) = [to - E(k) + iesgniu)]- 1 = [u - e{k) + itsgn(u)}- 1 
-oo 

(26) 

where Co = u — // and e(k) = E{k) — \i is the energy measured with respect the chemical 
potential \i > 0. 

a-The Averaged Single Particle Green's Function 

Using the free Green's functions G^(k,u) we compute the self energy as a function of 
the statistical average of the disorder potential. 

Z(k,u) = v 2 J ^(([/ sc (fc-^^-fc)))f||(fc,^f||(^fc)G(°)(p ) a;) 
f d 2 p 1 

= D sc v 2 J J^^ cos2 l2^ X ^ ~ *(p)][w - e($ +iesgn(u)]~ 1 , 

(27) 

10 



We replace the integration with respect \ by a polar angle integration on the Fermi surface 
E(k) = p. 

d 2 p _ pdpd8_ _ (ii+e)de / dx(p \-l\ ^ dx(p) ^ (n+e)de d\{p) 
(2tt) 2 2tt 2tt 2tto 2 V dd > \p=k F 2n ~ 2tto 2 2vr 

We obtain for the self energy: 
Z(k,cu) = D sc £ hA " ( ^J 2 dt J ^cos[\{ X {k)-x{k))][uJ-e{p)+iesgn{u)r l (28) 

E(fc, u;) = Sfl(fc, w) + zE/(£, a;), 

a,) = + M log( ^ A " ^ A 

47T TVJJ + jl 



— [i = vh\p\ = e(k); [i = vfikp 



(29) 



The real part of the self energy allows to compute the wave function renormalization Z. The 
value Z is determined by the band width cut-off A, the chemical potential ji and the elastic 
mean free path /: 



r 1 1 1 1 - ! 

Z=[l-0al!fl(w,/i)|a = o] = [1 - ^^t 1 + vha _ 1 ~ Lo 9(^a^I^ 

F M A 4 



(30) 



Considering the experimental situation we have — « 3 and M > 1 the wave function 



iiciitai aiiuauun we iiavc 

renormalization Z « 1. 

The averaged Green's function G(p, w) is expressed in terms of the renormalized wave 
function Z, renormalized velocity i>.r and renormalized life time Tr. We find: 

w) = 7-r, e*(p) = /mTr • p, w fi = uZ, tr = rZ" 1 . (31) 

u-t R {p) + i{^)sgn(uj) 

Here G +) ij(efl(p),a;) = — — f= .z_ is the retarted Green's function and G-^e^p),^) = 

^- — y- is the advanced Green's function. 

b-The ladder diagram-renormalized charge and transport scattering time 

Next we compute the effect of the ladder diagrams. The ladder diagrams renormalized 
the charge e and causes the replacement of the elastic scattering time r with the transport 
time r tr \ The calculation will be done at zero momentum. We introduce the renormalized 
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charge e^(fi) at a finite frequency fl. The current is determined by the product ev and the 
photon vertex W\(k, k). In our case the velocity normalizes separately such that the product 
vC^(k)C(k) = VRC^ R {k)CR{k) is invariant. The renormalized charge satisfies the following 
integral equation. 

^M/ 1 (p,^[((^ c (fc-^c/ sc r-fc)}}] W 2 f||Ap)f||(p,fc) 

(32) 



In order to solve the integral equation we use the explicit form of the scattering vertex 
product, f \\(k,p)T || (p, k) = cos 2 [|(x(fc) — x{p))]- 111 addition we replace the photon vertex 
Wi(p,p) = cos[x(p)] by cos[x(p)] = cos[x(k) + {x{p) ~ x(k))], using trigonometric relations 
we find: 

/ ^W 1 {p,p^t ll (k,00,k) = W 1 (k,k) f ^cos[xb1-p)]cos 2 [i(xW-xb-))] = ^i(^fc) 

(33) 

As a result we factor out Wi(k, k) on both side of the integral equation. We perform the 
disorder average [{(U sc (k — p)U sc (*— k)))] = D sc , introduce the notation for 1< ^ l ^ sc an d and 
perform the energy e(p) integration: 

Using this result we obtain the solution for the integral equation: 

e e 



i(n)D sc 1 _ i i 

1 4 2 l+if^rZ- 1 



e K (fi ->■ 0) ->■ 2e 



(35) 



The enhancement of the charge by a factor of 2 can be interpreted as increase of the transport 
time with respect the life time r tr ' = At. 

c-The contribution of the maximal crossed diagrams 

Following [16| we compute the maximal crossed diagram for the conduction electrons 
with a fixed chirality using the scattering vertex T\\(k,p). Following the discussion given in 
section III we obtain that the maximal crossed will become negative due to the identity 
x(— p) = x(p) + n - We introduce the notation for two particles Cooperon scattering: 

V[-p + Q,F,-k + Q,k-,n] = C[p,k-Q,Q] (36) 
12 



where —k + Q,k are the momenta of the incoming particles and — p + Q,p are the momenta 
of the outgoing particles. In the second notation Q represent the total momentum conserved 

— * — * — * 

in the process. We will work with the short notation C\p, k; Q,fl\. We replace [({U sc (k — 
p)U sc (p — k)))] by the statistical average D sc and obtain the integral equation: 

d 2 p 



C[k,k;Q,n] = D 



(2tt) = 



t\\(k, -p)t n(-p, —k + Q)f\\(-k + Q,p + Q)f „(p + Q, k)G + , R (-p,ou)G-, R (p + Q,cj + n) 
d 2 p 



+D S 



(27T) 



; C[k,p; Q, Q]G + , R (-p,u)G-, R (p + Q,uj + <>)!' (/.••//)!' i k + Q, -p) 



(37) 



C[k, k; Q, Q] is negative! In particular we observe that the n order term is proportional to: 



1 n-2 1 

~ j sin 2 [x(>) - x(Pi)} 11 cos2 [o(^(^) _ x(Pi+i)} 



(38) 



We solve iteratively the integral equation and perform the integration with respect the 
variables x(Pi)jX(P^)j--x(Pn-i)- We introduce the notation: 



I(fl,Q)D sc D sc /i f de(p) 



(39) 



where I = vt — v r t r . We find that C[k, k; Q, Q] is given by: 



C[k,k; 0,n] = -±D 8C [ 



+ 



W,Q)D 8 



= \d sc 



1 - 



i(n,Q)D s 



-] 



(40) 



This equation has the same form as we have shown in equation (24). Next we substitute 
the explicit result for J ("'Q) D " c given in eq.(39) and find: 



C[k,k;Q,n] = ^Dsc 



1 - 



-iVLrZ" 1 + \l 2 Q 2 



(41) 



This result shows that the Cooperon has a negative sign therefore we will have anti- 
localization. 



VI Computation of the conductivity 
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In this section we will compute the conductivity a(q = 0,Q — > 0) using the maximal 
crossed diagrams. We wind that the conductivity is determined by the longitudinal correla- 
tion function Ri t i(q, t — ti): 



-j.e[t-h](-ev) 2 



d 2 k 



d 2 p f d 2 ki 



d 2 pi 



(2tt) 2 J (2tt) 2 J (2vr) 2 J (2tt) 2 

w 1 {lp^w 1 (k 1 ,p 1 )5^k-p-^ 2 %-p 1 + 

(42) 

We perform the Fourier transform with respect to frequency Q and find that the conductivity 
is given by a(q = 0, Q) = Rl ^f- ' n ^ = a(Q). Using the averaged Green's function over the 



random potential we obtain: 



—i{ev) 



d 2 p f d 2 pi 



(2 



7T 



(27T) 



Wi(p,p)Wi(pi,pi) 



dse 



6[s]G(p,pi, s)G(p 1 ,p; -s) 



6[a){G&fc-a)G{fc,fia)y 



(43) 



Next we introduce G+ and G_, the retarded and advanced averaged Green's functions 15]. 
a^\Vt) represents the real part of the conductivity, computed with the averaged single 
particle Green's function G(p, u) and the ladder effect will be included with the help of 
the renormalized charge e R given in eq. (34). cr < - 1 - ) (fi) represents the correction to the 



conductivity obtained from the maximal crossed diagrams 

" emax de 

- I 

2tt 



16|. <j( )(f2) is given by: 



cr (0) (n) 



e R^R 



e 2 R ^R 
2h 



—i 



Z7T 



2tt 





— 2G +jR (e R , u)G- tR (e R , u + Q) 



-n 



4(0) 



2h 



l-i(&T)Z- 1 
1 + (Qr) 2 Z- 2 



(44) 



In order to compute the quantum corrections we need to compute the maximal crossed 
diagram. Using the maximal crossed diagram we consider the quantum correction to the 



conductivity fo 



B, 



lowing the so called Hikami boxes [7]. There are three such diagrams (see 



figure 4 in ref. 

(jW(O^O) 



The leading contribution to the conductivity is given by: 
i{ e R v ) 2 ^ f d 2 p f d 2 Q f° du. 



hiVt 



(2tt) 2 J (2tt) 



2n 



W 1 (p,p)W 1 (-p + Q,-p + Q) 
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- }R (p, u)G-, R (p, u)Gr +)R (-p + Q, a>)G_ H (-p + Q, u) vC[k, k;Q,Q = 0}v 



(45) 



Next we use the fact that the angular average in the limit Q = is given by (W\ (p, p)W\(—p+ 



Q, -p + Q)) 



Performing the momentum integration with the help of the residue 



theorem and neglecting the additive factor (see eq. (41) ) gives: 



e% f i dQ el , r 



<7 W(n->0) = ^ / ^r = ^Ln 
fin J J- Q hm 



I 



(46) 



This result shows tat anti-localization as obtained. The corrections depend on the phase- 
coherence length and the elastic scattering length. The phase-coherence length 
is limited by the length of the sample L, < L. The conductivity is given as a sum of 
the two parts,a(fi = 0) = a<®(Q = 0) + <7«(ft = 0). 



4e 2 

a(n -»• 0) « -^r-k F l 



1 + 



Tikpl 



-Ln 



L / JJ 



The real part of the conductivity as a function of the frequency ^ 



(47) 

Vtp = 4r), phase- 



coherence length and renormalized charge ^r(-^ 



is given by (34). 



2/i 



1 



l + ^M^- 1 ) 2 



4rck F r 



-Ln 



[£k F iz-y + \ 



^k F iz- i y 



+ 



If_L)2 



(48) 



In figure (1) we show the conductivity as a function of frequency for j- = 10 4 and Z » 1 
for different values of kpl = 1, 2.5, 5, 10, 20. Figure (1) shows that the conductivity increases 
with the decreases of kpl The figures also indicates that when j- — > and Q = the 
conductivity diverges. This results are in agreement with with the anti-localization theory 
predicted by the symplectic ensemble 



IV. Possible Experimental verification of the theory 



One of the difficulties in measuring the surface conductivity is the fact that the bulk 
electrons contribute to the conductivity. The bulk electrons form a three dimensional TI 
characterized by the inverted gap A. For the case that the chemical potential u > A 
the "' bulk" TI will form a conduction band. For this case the spinor will be given by 
*(r) = v^Ea^Et^fW^f^ In contrast to the surface case the particles carry 
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FIG. 1: The conductivity as a function of the frequency for three cases: kpl = 20 (the lower graph), 
kjrl = 10 kpl = 5 ,kpl = 2.5, kpl = 1 (the upper graph). 

two spin polarizations. As a result the Hamiltonian for this case will be replaced by: 

Hm = Y,^v(^wr^)-^)ci(k)c a (k)+Y J Y. E E vci(k)r a , a ,(k,p)CAP), 

<r=t4 % k p a=-t,i <T'=t4 

(49) 

where r CT)0 ./(fc, jo) is the scattering matrix. Using a model of parallel resistors we find that 
the conductivity will be dominated by the bulk contributions. In order to observe the 
surface conductivity we need to be in a situation where /x < A. For this case the bulk will 
be insulating and the conductivity will be dominated by the surface part. For this case 
electronic Raman scattering might be a good tool to observe the surface conductivity. This 
might work for laser frequencies ujp which excite the surface but not the bulk A > up > /x. 
In order to evaluate the Raman intensity we integrate out the anti-particles and obtain the 
effective two photon Hamiltonian: 



H (ext.2) 



d 2 k f d 2 q f d 2 q' f du { (ev) 2 W^ a \k + q,-k)W^ p) (-k,k - f) 



(2tt) 2 J (2tt) 2 J (2tt) 2 J 2tt 



C\k + q,u + u L + LO S )C(k- q,u)Ai(-q, -u L )Aj(-q, -lo s ), 
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(50) 



where Ai(— q, — ul) is the laser vector potential in the i direction, Aj(— <f, — u s ) is the 
scattering laser frequency and W i {k + q,—k) represent the photon vertex which cou- 
ples to electrons and holes. Taking the limit q — > we observe that the Raman fre- 
quency Q = ojl — Us can provide information about the electronic conductivity. The Ra- 



man spectrum measured in 19j] is given by: S(Q) = -[1 + (e^ n — l) _1 ]/m.D[fi] where 
D[tt) = -p- J °° dte 1 ^^^ ([p(t), p(0)]) represents the chiral surface electronic density- density 
correlation. Since Im.D[Q] is proportional to the conductivity <j(Q), Im.D[Q] oc Qa(Q), the 
Raman spectrum contains the information about the conductivity. For light the momentum 
q is negligible and we can replace (k + q, —k) ~ Ui(k), Wi a ' p \—k, k + q) ~ U*(k) which 

represents the electronic spinor projection in the i direction. Measuring the scattering light 
in the j direction for an incoming light in the i direction will be proportional to the density- 
density correlation (i and j correspond to the Cartesian coordinates on the surface of the 
TI, the incoming light vector k is almost perpendicular to the surface): 

D jti [qM = dte™ j ^ j ■^^U j (k)U i (k)U*(p^U*(p^(C\k+q, t)C(k, t^ip+q, t)C(p, 0))]. 

(51) 



The chiral properties of the surface electrons are probed with polarized light [19[. The 
circular polarized light e v , 77 = ± is expressed in terms of the linear polarization e 1 (k), t 2 {k). 
The chiral correlation function D_ t+ [q, Q] is given in terms of the linear polarized light and 
spinor s Ui(k): 

D_ i+ [q, Q] = {D XtX [q, Q] - D y Jq, Q]) + i{D Xty [q, fi] - D ytX [q, fl]), 
Ui(k) = bi. x sin^A^, k y ] - 5 i>y cos[x(k x , k y ] 
Im. Yl D iM^ °> fi ] K ^o-(fi)- 

i=x,y j=x,y 

(52) 

In Fig. 2 we plot the function D[q — > 0, f2] oc Qa(Q) using the conductivity given in Fig. 2. 
kpl = 15 (the lower graph), = 10, kpl = 5(upper graph). The line shape in Fig. 2 is in 
qualitative agreement with the Raman line shift reported in ref [3]. 

VIII Conclusions 



17 



Re[D(w/w max ,w max> =k / v,k / l)/(e 2 /(h *k / l) 




FIG. 2: The Raman shift £>[<f -> 0, fi] oc Oa(fi) for M = 5, M = 10 , k F l = 15 

To conclude we have computed the conductivity for the conduction electrons on the 
surface of the TI in the presence of a large chemical potential. Due to the TI topology 
the conduction band contains a half vortex at k = 0. As a result the return probabilty 
vanishes and the quantum corrections are given by anti-localization effect as predicted for the 
symplectic ensemble. Using this theory we discuss the implications to the Raman scattering 
recently observed in Bi 2 Se^. 
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